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We discuss the effect of matter reionization on the large-angular-scale anisotropy 
and polarization of the cosmic microwave background radiation (CMBR) in the stan- 
dard CDM model. We separate three cases in which the anisotropy is induced by 
pure scalar, pure tensor, and mixed metric perturbations respectively. It is found 
that, if reionization occurs early enough, the polarization can reach a detectable level 
of sequentially 6%, 9%, and 6.5% of the anisotropy. In general, a higher degree of 
polarization implies a dominant contribution from the tensor mode or reionization at 
high redshift. Since early reionization will suppress small-scale CMBR anisotropics 
and polarizations significantly, measuring the polarization on few degree scales can 
be a direct probe of the reionization history of the early universe. 



Subject headings: cosmology: cosmic microwave background — cosmology: observa- 
tions 



1. INTRODUCTION 

The large-scale anisotropy of the cosmic microwave background radiation (CMBR) mea- 
sured by the DMR onboard the Cosmic Background Explorer sateUite (Smoot et al. 1992) 
may be induced by density perturbations (scalar mode) and primordial gravitational waves 
(tensor mode) via the Sachs- Wolfe (SW) effect (Sachs & Wolfe 1967). If the anisotropy 
is mainly due to the scalar mode, the COBE result combined with observations of the 
large-scale structure would provide us with an important clue for discriminating between 
different cosmological models. However, it was recently argued that this anisotropy might 
be dominated by the tensor mode (Krauss & White 1992). Furthermore, it was shown that 
the tensor- mode dominance actually occurs in certain inflation models (Davis et al. 1992). 
Therefore, distangling the scalar from the tensor contributions is not only very important 
for testing the inflation model but is also needed to understand the formation of large-scale 
structure. It was then suggested that by comparing large- and small-scale anisotropy mea- 
surements, one can separate the scalar- and tensor-mode contributions (Davis et al. 1992; 
Crittenden et al. 1993). 

However, early reionization of the universe may influence the formation of structures 
(Couchman & Rees 1986) as well as damp out small-scale CMBR anisotropics (Vittorio & 
Silk 1984; Bond & Efstathiou 1984). The Gunn-Pctcrson test indicates that the universe 
must have been highly reionized by a rcdshift of five or greater (Gunn & Peterson 1965). 
Recently, the likelihood of early reionization by radiation emitted from young galaxies has 
been investigated in detail. In a standard cold dark matter (CDM) model, typical parameter 
values predict that reionization occurs at a redshift around 50 (Tegmark, Silk, & Blanchard 
1994); similar results were also obtained by numerical simulation (Fukugita & Kawasaki 
1994). When normalized to COBE observations, the CDM model very likely has reionization 
at redshifts around 28-69 (Liddle & Lyth 1994). With reionization at redshift ~ 50, the 
Doppler peak on degree scales may be damped almost completely away while large-scale 
anisotropics remain reasonably unaffected (Sugiyama, Silk, & Vittorio 1993). Indeed, it 
was attempted to have early reionization to smooth out excessive temperature fluctuations 
on degree scales predicted in the CDM model, in order to reconcile the model with the 
lowest hmits from South- Pole 91 data (Gaier et al. 1992). However, recent small-scale 
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anisotropy measurements from different groups have already hinted that the Doppler peak 
seems to be present at around the correct height for models based on adiabatic density 
perturbations without reionization. For nonstandard CDM models, some of them could have 
reionization occuring at rcdshifts early enough (Liddle & Lyth 1994) to suppress degree-scale 
anisotropics (Kamionkowski, Spergel, & Sugiyama 1993; Tegmark & Silk 1994a). Therefore, 
future anisotropy measurements on degree and subdegree angular scales would be crucial 
for determining the reionization history of the universe. 

Polarization of the CMBR is another clue that could have a great potential of probing the 
reionization history of the universe. Anisotropic radiation acquires linear polarization when 
it is scattered with free electrons (Rees 1968). There have been several works on calculating 
the small-angular-scale (< 1°) r.m.s. polarization of the CMBR induced by adiabatic density 
perturbations in an universe with standard recombination (Kaiser 1983; Bond & Efstathiou 
1984), and in a reionized universe (Bond & Efstathiou 1984; Nasel'skii & Polnarev 1987). It 
was shown that roughly 10 — 20% of the CMBR anisotropy is polarized on arc-minute scales. 
A thorough numerical calculation for both large- and small- angular-scale polarization of the 
CMBR induced by adiabatic and isocurvature density perturbations with standard recom- 
bination has been performed (Bond & Efstathiou 1987). Their calculations have confirmed 
earlier small-scale results and shown that large-scale polarization is insignificant. 

An analytic estimation of the quadrupole polarization induced by scalar and tensor metric 
perturbations was made in various cosmological models including matter reionization (Ng 
& Ng 1993). The r.m.s. temperature anisotropy and polarization of the CMBR induced by 
the tensor mode perturbation of arbitrary wavelength was also computed (Polnarev 1985; 
Frewin, Polnarev, & Coles 1993; Ng & Ng 1994). It was shown that large-scale polarization 
is greatly enhanced by early reionization. 

Similarly, the large-scale polarization of the CMBR induced by scalar and tensor modes 
within inflationary models was investigated under the assumption of no reionization (Harari 
& Zaldarriaga 1993). A detailed numerical calculation of scalar and tensor contributions 
to the CMBR polarization power spectrum in inflationary models has also been carried out 
(Crittenden, Davis, & Steinhardt 1993). Their calculations show that the polarization can 
reach a 10% level of the anisotropy in an universe with no hydrogen recombination, and may 
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be a useful discriminant for determining the ionization history of the universe. 

In this paper, we will investigate in detail the CMBR polarization induced by scalar 
and tensor modes in the presence of early reionization at redshift around 30 ~ 90. In this 
scenario, CMBR anisotropics on degree and subdegree scales are suppressed significantly. 
Hence, although the degree of polarization is still about 10%, the absolute polarization would 
be suppressed. In contrary, large-scale anisotropics measured by COBE remain reasonably 
unaffected while the degree of large-scale polarization is greatly enhanced. Thus, measuring 
large-scale CMBR polarization would become more important. 

In our previous paper (Ng & Ng 1994), we have given a detailed numerical calculation of 
the polarization of CMBR induced by pure tensor modes in an universe with and without 
matter reionization. It was shown that future polarization measurements with windows for 
the power spectrum from Z = 2 to Z = 50, at sensitivity of 10% level of the anisotropy, might 
have a chance to detect the CMBR polarization if early reionization took place at redshifts > 
90. Here we will combine these tensor- mode results with the scalar-mode contribution. The 
logic of this paper is more or less similar to the work by Crittenden et al. (Crittenden, Davis, 
& Steinhardt 1993). The difference is that we will not pursue models with no recombination, 
since they are disfavored (Tegmark & Silk 1994b) by the COBE FIRAS data (Mather et al. 
1994). 

2. METHODOLOGY 

We shall use the units c — h — 1 throughout. Our calculations are based on the standard 
CDM model with a flat metric: ds'^ — 0^(77) (^77^ — rfx^), where a{rj) and drj — dt/a{t) are the 
scale factor and conformal time respectively. Here we normalize the conformal time to unity 
today. In this metric, f^totai = ^cdm + = 1) where Qcdm and Qb denote respectively 
the cold dark and baryonic matter. The Hubble constant is Hq = 100/i km s^^ Mpc~^ with 
h — 0.5, and the baryon density according to nucleosynthesis is given hy flsh^ — 0.0125. We 
will approximate the ionization history by a step function: the universe is completely ionized 
before hydrogen recombination at redshift z^ec — 1375, so is after early reionization occuring 
at redshift Zi. Between is an universe with no optical opacity. This sudden approximation 
for the recombination and reionization processes is sufficient for the present consideration as 
long as their widths in redshift are reasonably narrow. The width will affect only calculations 
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on small-scale anisotropies. Since we concentrate on large-scale (> 1°) effect of scalar and 
tensor perturbations on CMBR, we include only the dominant SW effect (Sachs & Wolfe 
1967). 

To study how polarized photons propagate in the expanding universe, one need to solve 
the equation of transfer for photons (Chandrasekhar 1960). In general, arbitrarily polarized 
photons are characterized by four Stokes parameters, n = (n;, n^, n„, riy), where n — rii + rir 
is the distribution function for photons with I and r denoting two directions at right angle 
to each other. The equation of transfer for an arbitrarily polarized photon is governed by 
the coUisional Boltzmann equation. 



' d d \ 1 dn dhij 

dr] dx. J 2 dlnu dr] 



(1) 



where ut is the Thomson scattering cross section, N^. is the number of free electrons per 
unit volume, {]i — cos9, (p) are the polar angles of the propagation direction e of the photon 
with a comoving frequency u, and P is the phase matrix for Thomson scattering. 

The first term on the right-hand-side of equation () is the SW effect. For scalar modes, 
we have hij — J dkSkC^^'^kikj, where k is the wave vector and Sk is the fluctuation amplitude 
which satisfies the time evolution equation for density perturbations. Since we are consider- 
ing large-scale effect, we take 5k oc k'^rj'^, which is a good approximation for long-wavelength 
modes. As to tensor modes, we have hfj = J dkhke'^'^efj, where efj denote the gravitational 
wave polarization tensors, efj — e^ej — e*e*, and efj — eje| -|- e*ej; where and e* are two 
mutually orthogonal unit vectors perpendicular to the wave vector k. Then, the gravita- 
tional wave amplitude hk is governed by the equation of motion: hk + 2ahk/a -\- k^hk — 0, 
where dot means d/drj. The scalar and tensor power spectra are P'^^^k) oc T{k)k'^''~^ and 
P'^'^\k) oc fc*^*"^ respectively, where the power indices = 1 and = correspond to strict 
scale invariance. For cold dark matter the scalar-mode transfer function is (Bardeen et al. 
1986) 

y^^^ ^ [ln(l + 0.146fcr/,,)/(0.146%,)]' ^3) 

[1 + 0.242A;r/e, + {kr]eqf + (0.340A;r/eJ3 + (0.417A;r/e,)4]^ ' 

where r]eq is the time at which the energy density of radiation is equal to that of matter. 



The solution n for the equation of transfer assumes the form n = no + no(5n/2, where 
no and 6n are the unperturbed solution and perturbation respectively. We expand 6n = 
J dk n'e*'''^, where n' = aa + (3h. For the scalar-mode solution, the Stokes components Uu 
and riy both decouple from ni and n^, and it suffices to consider only the first two components 
of n with a = (1, 1) and b = (1,-1). Substituting the solution n and the Fourier expansion 
for hij into equation (), and expanding a and (3 in terms of Legendre polynomials, 

I 

/?(/x) = ^(2/ + l)APK/"), (3) 
I 

we obtain two coupled differential equations for a and /3, 

a + ikfia = ^T^k)k'^7]{l + 2P2) - q[a -a^- ^P2(«2 - /3o + /32)], (4) 
$ + ikfiP = -q[P - 1(1 - P2)(/3o - «2 - P2)], (5) 

where q = aTN^a. These equations are equivalent to equations (la) and (lb) of Bond & 
Efstathiou 1984 after making the approximations that the metric perturbations h = h^^, V 
and the baryon velocity v = 0. These approximations are valid as long as large-angular-scale 
calculation is concerned. 

Equations (^) and (^ can be casted into a system of coupled differential equations, 

do = —ikai + -T'2{k)k^r] 
di = —-k{aQ + 2«2) — qai 

d2 = ~k{2a, + 3as) + ^T"^{k)k''r] - ^(9^2 + /3o - (i^) 
5 15 lU 

/5o = -ikl3i - |(/5o + /32 + a2) 
/9i = -^A;(/3o + 2/32)-gA 

/32 = -^fc(2A + 3/33) - ^(9/^2 + l3o- «2) 
5 lU 

for I > 3, 

% k 

ai = -qai - _^ ^ + (/ + l)ai+i] 

A = -gA - m-i + + 1) A+i] • (6) 
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The coefficients ai and (3i are tfien solved by numerical method. 

In case of tensor mode, only decouples, and thus we choose the basis (Polnarev 1985): 

a = i(l - ^2)cos20(l, 1, 0) and b = + fj.'^)cos2(l), -(1 + /i2)cos20, 4^sin20), for the + 
mode solution. The x mode solution is given by the same expressions with cos20 and sin20 
interchanged. Defining ^ — a + (5, we obtain a system of coupled differential equations in a 
similar way, 

Co = -gCo - iHi + hk 

3 5 3 113 

$0 = -^0^^° ~ ^^^^ ^ ^^7^^ ^ 35'^^ ^ 10^° ^ 7^^ ~ 70^^'' 
for / > 1, 

ii ^ -q^i - ^^[i^i-i + {I + m+i] 

A = -qPi - + (/ + , (7) 

where the source term hk is obtained by solving numerically the equation of motion for hk- 
To describe the degrees of anisotropy and polarization, we compute the power spectra 
for the anisotropy, Cf, and polarization, Cf . To evaluate these functions, we expand the 
photon fluctuation distribution function in terms of spherical harmonic functions, i.e., 

Sn^J2 ^imYim, Siim = j 5vLYi*^(m. (8) 

l,m 

The total power spectrum is then given by (X^^ aj^a^^) = Ci' + C^, where () denotes the av- 
erage over all observation positions in the universe. For the scalar mode, it is straightforward 
to obtain (Bond & Efstathiou 1987) 

Cf^^^ = 87r(2/ + 1) I dkp(^)('^)l«f 

Cf^^^ = 87r(2i + 1) 1 cikp(^)(A;)|/3P|2. (9) 

The anisotropy and polarization power spectra for the tensor mode, C^^^^ and Cf*^^^ respec- 
tively, can be deduced in a similar way (Crittenden et al. 1993; Ng & Ng 1994). In deriving 
cf and C^^ , we have neglected the rotation from the /c-dependent basis to a fixed (lab- 
oratory) basis (see Bond & Efstathiou 1987 for details about the basis rotation). We have 
solved our transfer equations relative to this fc-basis and summed up all mode contributions 
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to the polarization power spectra in the same basis. Rigorously, this summation would not 
make any sense unless one has rotated the Stokes parameters from the /c-basis to the fixed 
basis. However, for Gaussian random perturbations, on the average, the expressions for 
Ci and Ci could well represent the degree of polarization in the laboratory (Bond & 
Efstathiou 1987). Having the power spectra, we can construct the correlation function, 

= ^^Ecr^W^^^/lcosG), (10) 

where Wi is the window function for detector, and is the separation angle. In actual 
observations the lower end of / is excluded by limited sky coverage, whearas the high-/ cutoff 
is fixed by the finite beam width. 

3. NUMERICAL RESULTS 

In this section, we present the scalar and tensor contributions to the CMBR anisotropy 
and polarization on large angular scales. We separate three cases in which the anisotropy 
is induced by pure scalar, pure tensor, and mixed metric perturbations respectively. To 
ensure that we capture the dominant contributions in the k integration when computing the 
multipoles in equation (), we have investigated the k dependence of certain Zth multipole. 
Roughly speaking, the dominant contribution to the Zth multipole comes from modes with 
k ~ / (note that r]^ = I). It is found that for / < 100, it is sufficient to consider the 
contributions from modes with k < 21, except for the tensor-induced polarization power 
spectrum. In the scalar-mode calculation, the k integration must be cut off at the point 
where linear perturbation theory breaks down. We thus set the cut-off at krec, where krec 
is the wavenumber which enters the horizon during the recombination era, since shorter- 
wavelength fluctuations probably have gone nonhnear (Abbott & Wise 1984). We refer the 
interested reader to Appendix for further discussion. 

Figure 1 shows the normalized CMBR anisotropy power spectra due to scale-invariant 
scalar (n^ = 1) and tensor {rit — 0) mode perturbations respectively. The constant behavior 
of 1(1 + 1)C^^^^'^'^^ /{2l + 1) for small I, with a scale-invariant spectrum in an universe with 
standard recombination {zi = 0), is evident from the figure. For the scalar mode, we have 
compared our numerical calculation for Zi = with another alternative approach which 
makes use of the SW integral formula (Sachs & Wolfe 1967). Using this formula, one can 
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easily derive an analytic free-streaming solution for the anisotropy power spectrum for Z > 
(Abbott & Wise 1984), 

dk 



cr«oc(2i+i)/|: 



+ :^krioSl - i''ji[k{r]o - rjrec)] 



I / + 1 ^ ^ 

jl-l[k{V0 - Vrec)] - TT^— rJZ+l[^(?70 - Vrec)] 



(11) 



where rjrec a-nd rjo are respectively the recombination time and the present time. Note 
that we have used rjrec ^ the lower integration limit in the SW formula instead of the 
usual decoupling time rjde- Under our assumption of the ionization history, they are in fact 
identical. We have plotted the power spectrum () represented by a long dashed curve in the 
figure. We see that, for small the result agrees very well with the numerical result. The last 
term in equation () is essentially due to a Doppler-shift correction for the world velocity of 
the source at the last scattering surface (Sachs & Wolfe 1967). In the numerical calculation, 
Thomson scatterings during the ionized stage before photon decoupling can effectively damp 
out this term. This effect is prominent mostly on small angular scales. That is why the 
high-Z C"^^^^s in the numerical result are lower than that in equation (). For the tensor 
case, similarly, we have plotted the power spectrum for Zi — caused by pure SW effect 
by putting the numerical solution of the evolution equation for gravitational waves in the 
SW formula (Ng & Speliotopoulos 1994), and found good agreement. Besides, late-time 
reionization {zi = 30 — 90) reduces the CMBR temperature anisotropy on small-angular 
scales, and shift the l{l + l)Cf^^'^^^ /{2l + 1) curves from scale-invariance. 

Figure 2 shows the normahzed CMBR anisotropy power spectra due to tilted scalar 
{rig — 0.85) and tensor [rit — —0.15) spectra respectively. It is evident that the constancy 
of l{l + l)C"^'^-''*-^V(2^ + 1) for small I is broken, however, the relative behavior of curves 
with different Zi values is similar to the scale-invariant case. Furthermore, the normalized 
magnitude of the power spectrum for large I is smaller when compared to the scale-invariant 
case. This can be easily explained by referring to the k dependence of the power spectrum 
in equation (). 

Figure 3 shows the polarization multipole to anisotropy multipole ratio for scalar per- 
turbations, Cj^^^^ /C^^^\ as a function of / with = 1 and 0.85. The ratio increases 
significantly for an universe which underwent an early reionization phase. In general, the 
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earlier the reionization takes place, the larger is the polarization-to-anisotropy ratio. Note 
that the ratio with early reionization has a peak around I ~ 10 — 40, which corresponds to a 
few degree angular scales. This makes the search for large-angular-scale polarizations more 
interesting. In particular, the peak of the Zi = 90 curve at / ~ 20 corresponds to about 
10% polarization in the CMBR fluctuations. Our results indicate that the ratio is rather 
insensitive to the variation in Ug- In fact, there is almost no discernible difference between 
the Ug — 1 and Ug — 0.85 results. To estimate the r.m.s. polarization-to-anisotropy ratio, 
[C''(0)/C"(0)]^/^, we sum / from 5 up to 50 with W/ = 1 in equation (), for spectra with 
Ug = 1, 0.85 and different reionization redshifts. This should correspond to typical large-scale 
CMBR polarization measurements. The ratios are listed in Table 1. In an universe with 
standard recombination, the polarization is much less than 1% of the anisotropy. However, 
the degree of polarization with early reionization at redshift ~ 90 is enhanced to 6.1%. 

Figure 4 shows the polarization multipole to anisotropy multipole ratio for tensor per- 
turbations, Ci^"^^ / C"^'^\ as a function of / with = and —0.15. The curves are similar 
to the scalar case, and the ratio is again insensitive to the variation in rij. Table 2 lists the 
r.m.s. large-scale polarization-to-anisotropy ratio. The polarization with early reionization 
at redshift ~ 90 for a scale-invariant spectrum is 9% of the anisotropy. 

So far, we have considered only pure scalar- or tensor-mode contributions to CMBR 
fluctuations. In general, these fluctuations can be induced by both scalar and tensor modes. 
In fact, inflation can generate both types of metric perturbations. In inflationary models, 
Ug is related to rithy Ug — rit + l, and there is a nearly model- independent relation between 
their induced anisotropy quadrupole moment: €2^'^^ /C2^^^ — —Int. For Ug — 0.85 and 
nt = -0.15, Cf'''^ ~ C^^'\ 

In Figure 5, we have plotted the total (scalar plus tensor) polarization multipole to 
anisotropy multipole ratio, cf^^'''^^'^'' /C^^^^~^^'^\ versus I for Ug — 0.85 and rit — —0.15, 
assummg C^^-^^ = C^^^l Also, in the flgure, we use short-dashed and long-dashed lines to 
represent the scalar and tensor portions respectively. The tensor contribution is dominant 
for almost all / in the standard recombination model, however, the dominance is taken by 
the scalar mode as Zi increases. Table 3 lists the r.m.s. large-scale total polarization-to- 
anisotropy ratio. The polarization with early reionization at redshift ~ 90 for this mixed 
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model is 6.5% of the anisotropy. 



4. COMPARISONS 

When comparing our result for the pure scalar mode with 2;j = in Figure 3 to Figure 
7 of Bond & Efstathiou 1987 and Figure 4 of Crittenden, Davis & Steinhardt 1993, we 
find that our ratio is about two orders of magnitude lower than theirs for small / and an 
order of magnitude lower for large 1. This may be due to the vanishing residual ionization 
subsequent to the hydrogen recombination that we have assumed in our whole calculation. 
In fact, we found that a residual ionization could raise the quadrupole polarization by two 
orders of magnitude (Ng & Ng 1993). This also explains why our r.m.s. polarization-to- 
anisotropy ratio (0.06%, see Table 1) for the standard recombination model with = 0.85 
pure scalar modes is much less than Crittenden et al.'s 0.4%. These discrepancies should 
be removed if we consider a more accurate model for the standard thermal history of the 
universe. Since the degree of polarization in the standard recombinaton model is already 
less than 1% which is well below the present detectable level, we will not pursue along this 
line and only concentrate on the reionization model. The 6.1% result with Zi = 90 in Table 1 
is, however, consisitent with their 7.9% result for the no recombination model. 

In Figure 4, the high-Z part of the curve with Zi — for the pure tensor mode is again an 
order of magnitude lower than the result in Figure 4 of Crittenden, Davis & Steinhardt 1993. 
But, the low-/ part is four orders of magnitude lower. Once again, the residual ionization 
can account for two orders of magnitude difference in the low-/ part. For the other two 
orders of magnitude difference, we suspect that it may be due to different definitions for 
the polarization power spectrum. Unfortunately, we cannot make an explicit comparison, 
since their paper did not show explicitly how to calculate the polarization spectrum. We 
sec from Figure 4 that the polarization to anisotropy ratio for / = 2 is 6.3 x 10^^. In fact, 
this value is close to our earlier analytic calculation for the tensor mode contribution to the 
polarization quadrupole moment in the instantaneous recombination model with vanishing 
residual ionization, which is about 9 x 10~^ (Ng & Ng 1993). 

As to the mixed-mode results, we compare our Figure 5 with Figures 1 & 3 of Crittenden, 
Davis & Steinhardt 1993. Once again, the solid curve with Zi = in Figure 5 is well below 
their result probably due to the reasons that we have mentioned above. The shape of the 



11 



solid curve with Zi — 90 agrees quite well with their no-recombination result. In particular, 
our value for the ratio of the total polarization quadrupole moment to the total anisotropy 
quadrupole moment (~ 5 • 10~^) agrees fairly well with their result (~ 8 • 10~^). Our r.m.s. 
total polarization-to-anisotropy ratio for the standard recombination is 0.13% (see Table 3) 
whereas theirs is 0.5%. With reionization at redshift equal to 90, we have 6.5%, which is 
comparable to their 7.4% for the no-recombination model. Although the general trend that 
the tensor contribution becomes subdominant as Zi increases is the same in both works, the 
details on how the dominance taken by the scalar contribution are different. In our work, we 
observe that when Zi = 0, there exists a window for / where the tensor mode is dominant. 
As Zi increases, the width of this window becomes narrow and disappears for Zi > 90 at 
Z ~ 5 - 10. 

5. CONCLUSIONS 

We have considered the effects of matter reionization on the large-angular-scale 
anisotropy and polarization of the cosmic microwave background radiation (CMBR) induced 
by scalar and tensor metric perturbations, with scale-invariant or tilted spectra. The results 
are rather insensitive to the power indices of fluctuation spectra. We separate three cases in 
which the anisotropy is induced by pure scalar, pure tensor, and mixed modes respectively. 
It is found that the polarization is insignificant in the standard recombination model. But, if 
reionization occurs early enough, the polarization can reach sequentially 6%, 9%, and 6.5% 
of the anisotropy. In general, a higher degree of polarization implies a dominant contribution 
from the tensor mode or reionization at high redshift. A 1% level sensitivity measurement 
could set constraint on the reionization redshift value, as well as provide information on 
whether the metric perturbation consists of a tensor component in models with a high red- 
shift reionization value. For instance, if polarization is detected to be 5%, then reionization 
occurs at a redshift between 45 and 90, rather independent of the types of the metric per- 
turbation and the power spectrum index. For the cases where Zi > 80, measurements with 
1% sensitivity can distinguish the scale-invariant tensor perturbation from the mixed mode 
metric perturbation. 

CMBR fluctuations have a 10% level polarization on angular scales less than 1° rather 
model-independently. In the CDM model with standard recombination, the predicted r.m.s. 
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I^T/T ~ 10~^ at 1°. Hence, to detect small-scale polarizations would require sensitivity at 
a level of AT/T ~ 10~^. In fact, this signal level can be achieved by using new technology 
and instrument design (Timbie, P. T. private communication). It is known that the universe 
has been reionized. Theoretical sides also predict an early reionization in dark matter mod- 
els. Our calculations show that, if reionization did occur at redshift ~ 90, the polarization 
to anisotropy ratio for either scalar, tensor, or mixed mode would have a peak of height 
of 10% around I ~ 20. This peak corresponds to an angle of about 9° {9 ~ tt/Z). When 
normalized to COBE/DMR anisotropy signals, this large-angle polarization is at a level of 
AT/T ~ 10^^ and would be detectable in the near future. Unfortunately, the early reion- 
ization would suppress small-scale anisotropics by an order of magnitude (Sugiyama, Silk, 
& Vittorio 1993), thus making small-scale polarization measurements elusive. It appears 
that measuring CMBR polarization on both large and small angular scales is a promising 
method for determining the ionization history of the universe. 

APPENDIX 

In this Appendix, we show how the /th multipole depends on the range of k integration. 
In Figures 6 and 7, we plot for the scalar mode the spectral power spectrum C"*-'^'' {k) , where 
Ci^^^ — J dk C^^^\k), as a function of k with Ug — 1, for Zi equal to and 90 respectively. 
In Figure 6, we also plot the dipole and quadrupole moments by using the SW formula 
(equation ()), and find very well agreements. The dipole moment (/ = 1) is sensitive to 
short-wavelength fluctuations. To ensure that we do not go beyond the linear regime, we 
set the cut-off k — krec in our calculation. We also investigate the k dependence of higher 
multipoles, and find that the dominant contributions arise from modes with k < 21 for the 
Ith multipole (however, it does not apply very well for small I). We illustrate this dependence 
for / equal to 2 and 20 in Figures 6 and 7. We see that reionization does not afi^ect much 
on the magnitude of the spectra but enrich structures at their high-A; tails. For / ~ 100, the 
moments will be damped out by reionization. Similarly, in Figure 8, we plot Cf^^\k) with 
Ug — 1, for I — 20, 40 and Zi — 0, 90. Obviously, reionization enhances the spectra by several 
orders of magnitude. The k < 21 rule applies very well for I > 20. Roughly speaking, the 
contributions drop by an order of magnitude when k is increased to 2k. 

In Figure 9, we plot C"^^\k) for the tensor mode as a function of k with rit = —0.15, 
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for Zi — and 90. We notice that the k < 21 rule apphes for high multipoles. The I — 20 
moment is damped by reionization. In Figure 10, we plot Ci^^\k) with rtt — —0.15, for 
Zi — Q and 90. Again, reionization greatly enhances the magnitude of the spectra. For an 
universe with no reionization, the spectrum falls off with k less abruptly. In this case, to 
capture the main contributions to the lih. multipole, one might need to include modes with 
k < lOZ in the calculation. However, this range of k is getting narrower for large I. At 
I ~ 100, it is sufficient to apply the k <2l rule again. 
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TABLES 



TABLE I. R.m.s. Polarization-to-Anisotropy Ratio for Scalar Mode 





[^{0)1! 
ic"(o)r 


iC"(0)J^ 




{Us = 1) 


{us = 0.85) 





0.0006 


0.00057 


30 


0.016 


0.016 


60 


0.040 


0.041 


90 


0.061 


0.061 



TABLE II. R.m.s. Polarization-to-Anisotropy Ratio for Tensor Mode 





Lc«(o)J 


^C«(0)J 




(nt = 0) 


{nt = -0.15) 





0.0025 


0.0019 


30 


0.022 


0.021 


60 


0.058 


0.053 


90 


0.090 


0.077 



TABLE III. R.m.s. Total Polarization-to-Anisotropy Ratio for Mixed Mode 









lc«(o)J 




{ns = 0.85, nt = -0.15) 





0.0013 


30 


0.018 


60 


0.044 


90 


0.065 
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Figure Captions 

Fig. 1 Normalized anisotropy power spectra. The solid and short-dashed curves corre- 
spond to the scale-invariant scalar (n^ = 1) and tensor {rit = 0) modes respectively, with 
reionization at Zi = 0, 30, 60, 90. The long-dashed curve is drawn by using the SW formula 
in equation (). The long-short-dashed curve is also drawn by using SW formula. For all 
curves in this figure and figures below, 0,b — 0.05 and h — 0.5. 

Fig. 2 Normalized anisotropy power spectra. The solid and dashed curves correspond to 
the tilted scalar (us = 0.85) and tensor {rit = —0.15) spectra respectively, with reionization 
at Zi = 0, 30, 60, 90. 

Fig. 3 Ratio of polarization multipole to anisotropy multipole as a function of I due to 
scalar mode perturbations. The solid and dashed curves correspond to the scale-invariant 
{us — 1) and tilted {ug — 0.85) cases respectively. 

Fig. 4 Ratio of polarization multipole to anisotropy multipole as a function of / due to 
tensor mode perturbations. The solid and dashed curves correspond to the scale-invariant 
{rit — 0) and tilted {rit — —0.15) cases respectively. 

Fig. 5 Total (scalar plus tensor) polarization multipole to anisotropy multipole ratio 
verus I, assuming equal scalar- and tensor- induced anisotropy quadrupole moments. The 
solid curves denote the total contribution. The short- and long-dashed curves correspond to 
the scalar and tensor portions respectively. 

Fig. 6 Spectral anisotropy power spectrum for scalar mode, C^^^\k), as a function of 
k with Ug — 1 and Zi — 0, for I — I, 2, 20. SW denotes the results obtained by using SW 
formula. 

Fig. 7 Spectral anisotropy power spectrum for scalar mode with = 1 and Zi = 90, for 
/ = 1, 2, 20. 

Fig. 8 Spectral polarization power spectrum for scalar mode C^^^\k) with = 1 and 
Zi = 0, 90, for I = 20, 40. 

Fig. 9 Spectral anisotropy power spectrum for tensor mode with rit — —0.15 and Zi — 

0, 90, for / = 2, 20. 

Fig. 10 Spectral polarization power spectrum for tensor mode with rit = —0.15 and 
Zi = 0, 90, for I = 2, 20. 
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